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abstract 

forces generated ^by^a e s?ngle en elLnd nV T h 831 '*^ ° f tHe aerod y namic 
simultaneous spinnine/whirlinT 8 ® 0 , g t ,abynnth seal executing a 
parameter mode? fo? a . nil* TT . been co ? d “cted. A lumped 
cavity with leakage is i\f i»i 8 T , SCa W *- tb coupling to an upstream 

technique. From this theory 6 i? °h 8 WUl L a " a PP ro P riate solution 
upstream cavity can in 1 Li Sh ° Wn that the P™*** of the 
direct damping by a factor of four a J*® mcnt the cross-stiffness and 
coupling are presented alone with he parameters that govern the 
simple uncoupled model i used ,S rf'T °" lhelr A 

for cross force generation Thk ° H d , fy lhe mechanisms responsible 
and lhe physical* significance of tl“ CCd W}™ is nondimensionalized 
Closed form algebra* fo™L/ r ' duCed PS™melers is discussed. 

cases. I. is als„ s h ow „ l L! ,he ',„,*', VC " f ° r f SOrae sira P" 
uncoupled model can hr rmrn. „ a lota , cros s-force predicted by the 

due to an inviscid flow with entrj swirl* a^T ° f 3 " * deal com P onent 
change in swirl crrmrH h r • ? • ■ d 3 v,SCOus Part due to the 

dependent ideal part is solely "resDonsTbl^f ^ >, 8,and ‘ Thc frequency 
damping. The facility designed IT? f ° r the rolord ynamic direct 

dependent forces s detr hed p U .‘° mCaSurc these frequency 

and usefulness of .£li/ Expenme " tal data infirm the validity 

calculating "h? d^mpin^ A ™*° d 

using only the stafir mroo. . ed on the force decomposition 

supporting the predicted cross Torre IS presen .‘ ed - Experimental results 
upstream coupling are' presented ^ 6 aUgmCntatlon due to the effect of 
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v = v*(i+n) 



pcfinition (units. SI) 

Direct damping. (N s/m) 

Cross damping. (N s/m) 

Non-dimensional direct damping. [1] 

Non-dimensional cross damping. [11 

Nondimensional sealing gap. [11 

Force acting on the rotor in the direction of 

instantaneous displacement. (N) 

Force acting on rotor in the direction tangent 

to the instantaneous displacement. (N) 
Frictional(viscous) cross force component. (N) 
Ideal cross force component. (N) 
Nondimensional depth of labyrinth knife. [11 

Depth of seal gland, (m) 

Depth of cavity upstream of seal, (m) 

Normalized energy carry-over factor.! 11 
Direct stiffness. (N/m) 

Cross stiffness. (N/m) 

Nondimensional direct stiffness. 11] 


Nondimensional cross stiffness. [1] 

Nondimensional labyrinth sealing pitch. [11 
Sealing pitch, (m) 

Axial gap from the vanes to seal inlet, (m) 
Static pressure in the labyrinth gland. (Pa) 
Total pressure upstream of first knife. (Pa) 
Static pressure downstream of seal. (Pa) 

Flow rate through first knifc/lcngth(kg/m s) 
Flow rate through second knife, (kg/m s) 

Gas constant for air. 287.15 (J/ kg K) 

Radial displacement of whirling shaft, (m) 
Magnitude of radial displacement, (m) 

Nondimensional spin rate. [1] 

Air temperature. (K) 

Time, (s) 

Swirl velocity inside the seal gland, (m/s) 
Nondimensional whirling frequency. [1] 


180 


0=62/8!* 

P = P 2 /Pi 
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Nondimensional divergence. [1] 

Kinetic energy carry-over coefficient. [1] 

Nondimensional swirl gradient parameter. [1] 

Nondimensional seal flow rate. [1] 

Sealing gaps for the sealing knives, (m) 
Nondimensional whirling eccentricity. [1] 
Darcy friction factor for the stator. [1] 

Darcy friction factor for the rotor. [1] 

Kinetic energy carry-over sensitivity. (1/m) 
Flow coefficients for the scaling knives. [1] 
Harmonic pressure perturbation. [1] 

Harmonic velocity perturbation. [1] 

Fluid density, (kg/m 3 ) 

Nondimensional swirl parameter. [1] 

Swirl parameter as observed in whirling 
reference frame. [1] 

Whirling frequency of the shaft. (1/s) 

Shaft rotational speed. (1/s) 

Steady solution for the case of centered rotor. 
Quantities at the first or second knife 
respectively. 


I. INTRODUCTION 

One of the major sources of rotordynamic instabilities in hiuh 

SnSSfUr ' ur h b0m !h hi " e 7 is asymmetric pressure dis'ribu.ions t 
seal 8 fL™ , f labyrmlh u sca,s which generate cross forces. If the net 

direction nf Tv J°, thC lnstantanc °us whirling displacement is in the 

direct' on of this whirl, instability will be promoted. 

l ^ possibility of labyrinth seals creating unstable rotor whirl has 
been known since the 1940's as reported by Den Hartog Ml Alford 121 

pr C d ic Ti 0 ° n ma 0 S f !; « "T lh ? «' sl “ P ro P 0 « analytic ,hi 

predict on of labyrinth seal forces. However, neither of these models is 

particularly useful, because the influence of the inlet swirl, which Is 
nown to dominate in the generation of cross forces, was neglected bv 

Ihe knives 8 * to ’S P ™ mct * T model that cou P Ie s the axial flow over 
inside the ceai one ' dimens, o nal continuity and momentum equations 
side the seal gland was developed by Koslyuk [4] The gland-demh 

assr „:; ,h d ir r r r d cccen,rici,y was 6l r 

whirl! no no .[ direct damping is predicted from this model for a rotor 
hiding with a parallel precession. Subsequently, Iwalsubo [5] added 

nc'uI'rSs'rdT j? T?"' , f °I ' hiS Var?a,ion Many o^c'r ’au ho , 

" „ a ", d , nouc m Fu j ikawa ’ *»">«*« aid Abe 

oj, acnarrcr and Clulds |9), and Martinez and Lee 1101 have used simil-.r 

lumped parameter models to predict the eross forcis for a gma" variety 
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of geometries and flow conditions. All of these models have assumed 
constant upstream and downstream boundary conditions. 

Analyses that consider variations within each seal gland, by 
using multi-control volumes within each seal gland, have been 
conducted [8,11]. However, this approach is cumbersome and the 
different control volumes must be coupled with ad-hoc assumptions or 
analysis. No advantages have been clearly shown from such models 
over the single control volume models. 

Computational Fluid Dynamics (CFD) methods have also been used, 
most notably by Nordmann and Wciser [12], to predict results similar to 
those obtained by the lumped parameter models. These codes may be 
very useful in modeling some of the detailed flow fields in the glands, 
which can be used for the sub-models to predict carry-over coefficient 
variation, friction factors, etc. However, it appears that the simple 
lumped parameter models contain the dominant fluid physics, which is 
asymmetric injection of swirl momentum, necessary to predict 
destabilizing forces, at least for multi-cavity seals. 

All of these analyses have assumed uniform upstream and 
downstream boundary conditions. These assumptions may not be 
adequate for the prediction of the rotordynamic forces in short seals, 
where the end conditions may greatly affect the perturbations in both 
pressure and the swirl component of velocity. 

The first experiments to measure the self-exciting forces in 

isolated labyrinths were those of Benckcrt and Wachtcr [13,14,15], who 
measured the static pressure distributions around the casing and 
integrated them to find the direct and cross forces due to a statically 
offset rotor. Many different geometries were tested over a wide range of 
flow conditions. They determined the general prcssurc/inlet kinetic 
energy scaling parameters for statically offset seals. Brown and Lcong 
[16], Thieleke and Stetter [17], Kanki and Morii [18] and Hisa, Sakakida 
and Asatu [19], have all made similar measurements yielding a good data 
base on the displacement dependent rotordynamic coefficients. 

Experimental investigations on the dynamic characteristics of 
labyrinths have been conducted by Wright [20], Kanemitsu and Ohsawa 
[21], Scharrer [22], Scharrer and Childs [23] and Millsaps and Martinez 
[24], 

Kostyuk-type lumped parameter models are capable of predicting 
the cross-stiffness coefficients for long labyrinth seals (more than 6 
chambers) to within 25% of the experimentally obtained values [10]. 
However, the situation with respect to short seals, and with the dynamic 
coefficients, especially the direct damping, has been far less 
satisfactory. The measured cross-stiffness coefficients in two and three 
gland seals, as measured by Benckert [12], are more than 100% larger 

than those predicted by theory in most cases. Measured dynamic 
coefficients for both long and short seals arc not well predicted by 
theory. Discrepancies between the measured and predicted direct 
damping of nearly 500% are shown by Scharrer [21] for some 
conditions. While the results from many computations have been 
reported using the various lumped parameter models, little physical 

understanding on the mechanisms that generate dynamic forces has 
been obtained from the models. In particular, no analysis has been 
presented that clearly delineates the importance of the various 
geometric and flow related parameters and explains, in physical terms, 
the origin of the damping forces. The general mechanisms that 
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1 trijz; 4 °z°r ls - ■ 

Which includes^ the cf^'^r^pstream 0 ,mr ? ducc a h . ncw analytical model 
responsible for the much h "32 P *’ *f ,Ch i§ bC ' ieVCd 10 be 
direct damping obtained cxner ° xpcctcd stiffness and 

mechanisms of the damping forces m labyrinth seafs'^ ^ " alUrC and 

uniform flow Tn the™ volume ^ m ° dd ' S prcscntcd ,ha ‘ allows for non- 
continuity and momentum cauafloi? 3 " 1 ° f ,hC whirlin S Coupled 
single gland seal with the added TeL H g ' VCn , f ° r bolh lhc wh *rling, 
for the non-whirling finhe volume ' ° ^ n0 "- Uniform inlet, and 

leakage path to a large center cavhv ThT ^ which has a 
flow significantly complicates the ™| T , bC add,tl °" of the extr a leakage 
introduces an e?,eJ^ this because it 

is presented, based on harmonic Hi' * approximate solution method 
The parameters that contml the ? gmg ’ l ° deal with lhis Acuity, 

identified. Results from thi " new Lr,®" , ' n “‘" 0n cross f °'«s are 

impact due to the non-uniform^ How 'S h "’ a ' S T a lar ^ 

with the seal perturbations. P f f,eId due t0 the coupling 

is discussed!* 1 °Thc mecSsms ’“£? ,t ^ f ° ,CCS in labyrimh sca 's 
damping will be delineated All the « cneral,on of rotordynamic 
developed with uniform inlei ennH .• th phySlt : a, arguments will be 
complexity. However all of ihe d ‘° nS l ° aV ° ,d unne ccssary algebraic 

readily generalized to account for thfTT ‘° be prCSCnlcd can bc 
reduced equations are nnnH* l ^ e ^P stream non-uniformities, These 
frequency-dcpcndcm di«c, “d ' a " d for lhc 

significance of the varil n, ra forces are given. The physical 

behavior provided The roirtrH C [? ,S d ‘ scussed > and the scaling 

be shown to be the sum of ordy " an ? lca,| y destabilizing cross force win 

an inviscid flow, and a "Wscous^ C ° mp0ncnls - A ? one due to 

damping will be 'shnwn • ° US par f due t0 frictional shear. The 

From this decomposition, a*' method* 11 ?* fr ° m ^ ' dcaI component. 

cncfficicntSs from purely static measurement arc ‘'showm **" dampinS 

methods used to dc.'crmmc 'these T*"' ■ ap ? aralus and measurement 

‘ns, rumen, a,i„„ o and ^^XSres ‘l™™. 

will be provided. AlthoS^rprecisl 116 the01 ? l ° thC experimenta I data 
to the experimental data is c ° mpar,son of lh e coupled model 

measurement of the axial eaD of the r dUC l ? the ,ack of contro1 or 
pressure perturbations to the Center ^huh^il Wh ' Ch t? S lhC upslrcam 
data support the use of the ideal/ hb p cn . um - Thc experimental 

ssrar z 

Bencher, for staticaMy offset ^ ° f 

2. MODEL FOR UPSTREAM COUPLING 

the assump'ion eri ,ha" la !he and bou a n n H! ! '* iCa ' r ? earCh P rnc “ d '« initially on 

labyrinth l *seaJ° were 'uniform ”* Thtis ' hc S la " d 

-y an upstream reservoir where the' press™' t? ^ « 
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spatially and temporally uniform, and discharges to a downstream 

volume at a constant static pressure. . , 

Indications that substantial upstream non-uniformity ma * 
and be of importance were provided from l ' m « rcs ® lvcd P - h 
measurements in the MIT Labyrinth Seal Test Facility (LSTF) wUh a 
whirling seal, just upstream of the first knife, and from steady pressure 
measurements \n the MIT Alford Force Test F acilUy CAF^F } upst. ream *nd 
downstream of a statically off-set seal, that was attached to a shrouded 

turbinC Th r c 0t ° r x imental values of cross stiffness and direct damping 
measured in the LSTF were consistently 2 to 4 times the values P^ dict *£ 
by the lumped parameter model. Several attempts were made 
reconcile the ? simple model with the experiments by parametencally 
varying the sub-model coefficients such as friction factors, dlsc | ia rgJ' 
coefficients, kinetic energy carry-over coefficients etc. It was f d 
that within reasonable ranges of variation, no combination oi the 
coefficients would yield results in agreement with the experimental 
data Therefore, the basic assumptions of the model were invest, gated 
This reassessment led to an extended model where t J f 
upstream flow, induced by the flow perturbations in the «!“ d '“ clf f 
included. This extended theory is capable of predicting the high 

levels that were measured. . . r , 

Before embarking upon algebraic manipulations, it is Evolved 
provide a qualitative description of the coupling mechanisms >nvolvcd 
When the shaft is offset and the scaling gap vanes around the seal 
ncrimeter a low pressure area will develop upstream of the seal in 
vicTity of the widest gap. Because of this, the velocity magnitude w.l 
be higher there, .ndV in particular, the all-important tangential 
velocity will have a relative maximum near the widest gap. When the 
fluid carrying this excess tangential momentum enters the seal cavity 
and mixcs'wi.h.he swirling seal flow, i, will preferential* energy « 
in this area. The result will be a positive pressure gradient QP/dQ), m 

this area, and hence a maximum P in the seal about 90 ahead of this 
maximum gap. This will then produce a forward whirling force For a 

concrete formulation of these effects, consider the geometry of ig 
1. This geometry reflects that of the LSTF and most other test section 
configurations that have been used to measure rotordynamic orces. 

The swirl vanes, which arc located ^upstream from the first knife, have 
an effective radial gap of 8y . There is some reduction of the effective 
flow area due the vane metal and boundary layer blockage. These vanes 

deliver air into the first cavity with an effective swirl angle of a v , 

which is the metal angle minus some small turning deviation. The 

cavity is h deep and is sealed from a large center volume by an axial 
face labyrinth seal with a gap of 5- . Since there is no net flow into this 

center cavity, the pressure here is uniform and the same as in the swirl 

cavity, namely P*. The one-dimensional continuity equation for the 

upstream cavity (swirl chamber) is 


+ r^ (Pil,h,Vi) + qi ' qv + qc>ou ‘ ' = ° 


( 1 ) 
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where \ is no longer constant, q v is the flow rate per unit length 

issu.ng from the swirl vanes and the q c 's are the radial flows in and out 

•siiffirilni f Cnt< ! r ? avi,y L ’ respectively. Incompressible relations are 

low These n Calm8 ‘ hcsc . nows sincc thc transfer velocities arc very 

low. These flows can be written as * 


t lc=Mc^i(P i -P*) 


( 2 ) 

This relation is fundamentally different from those for q v , q, and q 2 in 
ere is no ow to or from the center volume when thc seal is 

centered in the casing because P i= ^*. This introduces an essential 
nonlinearity into the analysis and hence must be dealt with in a 
erent manner. Likewise the momentum equation in this cavity is 

f(P.I i h l V i ) + .l.^(p i I j h i V i 2 ) + q 1 V 1 -q v V v 

+< lc.ou l V i - q c ,inX: + *.(21; + h;) — T f h; Q 

R * dQ ( 3 ) 

where V e is thc swirl velocity inside the center volume. In this cavity 
the cross sectional area, l.h,, and thc vane gap. C are constant. 
However he inlet swirl component of velocity. V v . is not. The angle of 
thc fluid leaving the vanes. a v . is constant. Therefore, a drop in the 
pressure at one location in this cavity will induce a greater mass influx 
and hence a higher swirl velocity at that location. This is tL essence of 
the mechanism that augments the forces. 

. ThC usual i' D continuity and momentum equations for the seal 

?ons1raTnts Pre of m the by 77 Au ' honi |5 ' 6 - 8 ' 9 - 101 ar <= still valid within thi 

additional term« e m ° dC V except that thcir linearization will yield 
additional terms from the upstream pressure and velocity non 
uniformities. The continuity equation is velocity non- 


afpi(h +s,)] i a 

dt + r, a© + )vj + q 2 ~qi =o 


( 4 ) 


and the momentum equation is 
9[pl(h +5,)VJ i a 

5 + T- a? ,pl(h + *' ,v! 1 + - 1. v . + '.' - '■(! » 2I» + = 0 (5) 

* R, o0 

and the q's, other than those to and from the center cavity, arc given by 






( 6 ) 


center cavity. However, as previously stated the nature 


185 


of the oscillating flow between the two upstream volumes requires some 
modifications and careful treatment. This is because these oscillating 
leakage flows vary as i/pi, where P- is the upstream pressure 

perturbation, rather than linearly with (first order) or P* (zeroth 
order) like the others. Thus, their effect can be disproportionally great 

at small offset amplitudes. 

First the zeroth order solutions are obtained for the centered 
rotor. The velocity, pressure and density in the upstream swirl cavity 

will be denoted by V*, P* and p‘ respectively and similarly in the seal 
gland, V*. P*and p*. The pressure and velocity in the swirl cavity are 
expressed as 

P, = P.-fl + ic 1 ' 8 - 00 ) Vj = V*(l + f|i C i(9_Q,) ) (7) 

and in the seal gland by 

P =P*(l + Se ,(e "° ,> ) V =V*(l + f|e i(e - Q,) ) (8) 

Where it is understood that the real part of all expressions are to be 
used. The perturbation expressions for £ it | and tj are substituted 

into the continuity and momentum equations for both the upstream 
swirl cavity and the seal gland. The nondimensional perturbation 
leakage flow into the center cavity is 

<9) 

The first harmonic component of this function will be extracted by 
averaging over one period. The first harmonic of the perturbation 
reduces, after some manipulation, to 

^Rc{| i e w }p=^^|cos I ^d4'=^^V2B^,|j (10) 

Where B is the beta function. From this the first harmonic of q cout “‘Ic.in 
is found to be 


% = 1.57377^ J-Ar*-4ir‘ /2 4 i e i(0 ' Q ‘ ) 


(ID 


Similarly the first harmonic of q Cinul V, ~qc.in V c is 


1. 57377 ?( Vc * + Vi ‘)^ ei(0 

M, l( p i- p 


a) 


(17) 
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Th f firs ! harm onic of the perturbation continuity equation for the 
swirl cavity, after dividing by q* is 


1 


L57377uX 

— i 

JV 

* 

to 

2 

^ 1 




4- 

k-if 

+ * 








(13) 


and the swirl cavity momentum equation after dividing by q*V* is 




^Prpi-HQVI* 

2y 


Ipilihi 





(14) 


The nondimcnsional continuity equation for the 
includes the effect of the upstream cavity coupling 


seal 

is 


gland, 


which 



and the coupled momentum equation for the seal gland is 


(15) 



A. r (l + 2h)(toR, - V*) 2 )+ 

.W U. ) R. 

1 + ^( X ’ 1 - M + 2HXV* - (OR.) 2 ) + p*ih^ - aji 


(16) 


C-3 
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This can be written compactly in matrix form as 


(7 

^ 1,2 

^ 2,1 

^ 2,2 

Z3.1 

^ 3,2 

^ 4.1 

^ 4,2 


^1.3 

7 \ 

^1.4 

^2,3 

^2,4 

^3,3 

^3.4 

Z4.3 

^4,4 j 


i, 


Rj 

Tli 


r 2 

k 


r 3 

T| 


R 4. 


(17) 


This is not a linear system. The elements Z M and Z 21 contain 1^1 1/2 . Of 

course if £j=f|j=0 the system decouples and the original 2x2 system for 
the isolated seal with no upstream coupling is recovered. To solve this 
system a simple iteration scheme was employed. Equation (17) is first 
solved with 8 C = 0 then from this solution, £{, is used on the left hand 
side to calculate an updated solution. This iteration procedure continues 
until a predetermined convergence criterion is satisfied 

l^ +, -|"l<e 08 ) 

Now this model shall be used for predicting the influence of this 
coupling mechanism. In addition to the parameters that control the 
generation of rolordynamic forces for an isolated seal, the parameters 
that characterize the influence of the upstream coupling arc 

V. The ratio of the swirl cavity area to the seal gland area, (ljhJ/Oh). 

2. The relative size of the axial sealing gap compared to the radial gap, 

5 c /s;. 

3. The swirl velocity inside the center cavity, V c . This is strongly 

influenced by rotation of the seal disk. 

4. The relative whirl eccentricity, r/8*. This is a purely non-linear 

effect. For the linear system all of the forces are directly 
proportional to the whirl eccentricity. 

The effects of the cross force augmentation of each parameter 

will now be considered separately. According to the model, if there is no 
leakage into the center cavity, the effect of the upstream coupling 

always acts to increase the magnitude of both the cross stiffness and 
direct damping, and in the same proportions. Figure 2 shows the ratio 

of the direct damping from the coupled model, with 8 C = 0, to that of the 

uncoupled one for various swirl chamber to seal area ratios. As the 

swirl cavity area approaches zero the predicted force augmentation does 
not vanish but, approaches a value of 1.62. This residual effect in the 

absence of the first cavity is due to the condition imposed at the swirl 

vanes. In the simple model V] is constant. If instead the vanes arc close 
coupled a reduction in the gland pressure will still bring in more flow 
and hence will induce a higher swirl component locally. The maximum 

increase in the cross stiffness and direct damping over the uncoupled 
model is about 4.42 and occurs at an area ratio, (lih,)/ (lh) of 1.35. Even at 
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an area ralio of 10 the forces arc increased by a factor of two. The force 
predicted by the coupled model asymptotically approaches the 
uncoupled one as (l,hj)/(lh) -» ». The coupled and uncoupled models match 
to within 1% for an area ralio of about 80. Well before this value the 
assumptions of the model probably break down. In particular, 
significant variations in the perturbation quantities arc likely to occur 
in the axial direction within the swirl cavity. 

The presence of the axial clearance between the swirl cavity and 

the large center volume permits for a "venting" that reduces the 

magnitude of £ it This effect tends to mitigate the large augmentation of 
the forces that the upstream cavity may induce. Figure 3 shows the 

direct and. cross force vs. the relative leakage area 8 c /8[. It is assumed 
that V c - V ; for simplicity. As 5 C goes from 0 to » both force components 
go from the fully coupled values to those predicted by the uncoupled 

model. However, this does not occur when V C *V*. The forces are very 
sensitive to small changes in the axial gap when it is less than 8* 

However, when 8 c /8, >1, there is a greatly reduced sensitivity to small 
changes in axial gap. 

The model predicts that the swirl velocity inside the center 

volume can have a minor impact, on the order of 10-20%, on the seal 
pressure perturbations. For cases where there is no seal rotation it is 
probably safe to assume that V c =0 . This is because the tangential 
momentum feed into the seal is of perturbation order and the shear 

stresses acting to retard the flow are of order unity. For cases with seal 

rotation, it would be very difficult to estimate the swirl velocity inside 

the center cavity. Figure 4 shows the effect that changes in the center 

cavity swirl velocity have on the forces. 

In the absence of the leakage flow non-linearity, the theory 
predicts that the forces should scale with whirl eccentricity and hence 
the rotordynamic coefficients, and should be independent of the whirl 
amplitude. The nonlinearity due to axial leakage is shown in Figure 5, 
which shows K„ a nd the relative eccentricity r/S*. The behavior of 
C,, is the same as for K ly . The direct force is much more sensitive to the 

whirl amplitude than is the cross force. At large whirl amplitudes, the 

predicted forces approach those obtained for 8 C =0 (i.c. the fully coupled 
case). However, as r->0the center leakage flow is able to "kill" the swirl 
cavity pressure perturbation completely. This effectively decouples the 

whirling seal from the upstream cavity. 


3. IDEAL - VISCOUS FORCE DECOMPOSITION 

It will now be shown that within the constraints of the lumped 
parameter model, the total cross force can be separated into two additive 
components. One of these is the ideal part, due to an inviscid inlet flow 
with entry swirl and the other is the viscous part, due to the change in 
swirl brought about by frictional shear stresses in the gland 

The following development shall be done without consideration 
of upstream coupling. However, all of the arguments can be readily 
generalized to account for such coupling. ^ 

„ nn , In . ord . e . r t0 isolale these two contributions it is useful to 
nondimensionalize the continuity and momentum equations for the seal 
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gland. Equations (15) and (16) can be written in terms of the following 
non-dimensional parameters as 




• \ 2 
0^2 

W>) 


oL 

yd 


(1-W) 




{"d~ = {(a _ + K + aL(1 " W)i } £ (19) 


and 




oLW L 
+ 


yD A 2 ap‘ 2 D 


1 + 


ctL f 
4DH ' 


X, - X r ^l + ~)(1 - S) j + ^(1 - W)ijf| = re (20) 


All of these non-dimensional parameters may be categorized as 
geometric, kinematic or flow related. The geometric ones arc 


e = 


A 

r 



5; 

sf 




(21) 


Note that K has been considered along with the geometric parameters 
because it is analogous to a convergent/divergent gap (Alford effect) in 
the effect that it has on the direct and cross forces. The kinematic 
parameters that specify the motion of the seal are 


S = 




( 22 ) 


and the dynamic parameters that indicate the axial flow rate or pressure 
gradient, the inlet swirl and the change in swirl respectively are 



r=i--4 (23) 

v 


There are two combinations of these parameters, or and a(l-W), which 
will be shown to be very important. 

Cramer's rule can be used to write an explicit expression for the 
nondimensional pressure perturbation which, upon integration, yields 
the forces acting on the rotor. The forces are proportional to P*£. In 
general the shear stresses may also contribute to the forces. However, 
from the equation solutions it has been shown that the shear forces due 
to velocity perturbations arc small when compared to the forces arising 
from pressure perturbations [24], For moderate whirl relative inlet 
swirl angles, which occurs providing the following condition is 

satisfied 


l + q 2 p 2 
A 2 


»— (1-W) 
yD 


(24) 
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the full expression may be simplified to yield the following real and 
imaginary parts for the pressure perturbation normalized by both the 
flow and eccentricity 


c K-fl-lV — (1-W) 2 

k V aJ P 

A £ " 2 o 2 L? 

1+o|i 


aL [ 


A 2 E 


r + (l-W)(-D-K + l 


-i)li 

a 


,+ap + 5v 


(25) 


(26) 


The first of these is real and therefore is in the direction of the 
minimum gap and hence is proportional to the direct force. The 
imaginary part is proportional to the destabilizing cross force. Both 
forces are proportional to P*A 2 . For low inlct-lo exit pressure ratios, this 
reduces approximately to a simple pressure difference scaling, Pi-P 0 , as 
one would expect for an incompressible flow. For this to hold, all of the 
other nondimcnsional parameters must be kept constant. In particular, 
the inlet swirl flow angle as observed from the whirling rotor must be 
fixed as the pressure difference is changed. As an example, if the axial 
pressure difference is increased the inlet swirl velocity and whirl speed 
must also be increased to maintain constant, 0 and W, if this scaling is to 
be used. This is a generalization of the relationship used by Benckert 

[12] for statically offset seals. 

The direct force is mainly due to the kinetic energy carry-over 

variations and differences in the nominal sealing clearances as seen in 
Equation (25). These effects generate a direct stiffness. The smaller 
terms, that arc proportional to a 2 (l-W) 2 , generate direct inertia 

coefficients. It is possible that this direct force may alter the natural 
frequency of the rotor slightly. Therefore, no further discussion of the 

direct force will be given due to the small impact it has on rotordynamic 
stability. 

The cross force can be seen from Equation (26) to be the 
summation of two terms, one proportional to o(l-W) and the other 
proportional to or. The nature of these two contributions and their 
crucial differences arc very important and will lead directly to the cross 
force decomposition to be given. The kinetic energy carry-over term, K, 
and the "Alford" [2] term due to seal convergence/divergence, (1-1/a), 
do not generate cross force in the absence of swirl. But do alter this 
force in the presence of swirl. A convergent seal and/or kinetic energy 
carry-over will tend to increase the cross force magnitude. This is 
destabilizing when W<1 or equivalently when the whirl frequency is 
less than V*/R,. However, this will enhance stability, in the forward 
direction, when Q is greater than V*/R,. 

The total cross force, Fp, will now be split into "ideal" and 
"viscous" contributions. The ideal or inviscid part, F 0 , is the cross force 
that would be generated in a purely inviscid fluid with inlet swirl. It is 
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proportional to cr(l-W), The viscous contribution, F r , which is 
proportional to al\ is due to the frictional forces changing the swirl 
velocity inside the seal gland. The two forces tend to be of the same 
magnitude but, the nature of these two contributions is quite different. 
The ideal component is dependent on the rotor whirl frequency and 
hence will contribute to the direct damping. The viscous component is 
independent of the whirl speed and hence can not contribute to the 
direct damping. The frequency dependent behavior of a seal in an 
inviscid fluid can be readily explained by considering a simple change 
of reference frame. Let the relative inlet swirl, denoted by V^Vj-ftR,, 
be the inlet swirl velocity as measured by an observer rotating in the 
whirling frame. The non-dimensional swirl parameter in the whirling 
frame will be similarly denoted by a and is 

o = a(l-W) = ^p- ( 27 ) 

The cross force generated by an inviscid flow at an inlet swirl of \ i at a 
frequency of ft is identically the same as the cross force that would be 
generated for a static offset with the associated relative inlet swirl of V i , 
This equivalence can be seen by noting that the governing continuity 

and momentum equations are invariant under Galilean transformation. 
As a particular case of this argument, consider a rotor whirling so that 

the minimum gap is traveling at the same speed as the inlet swirl. Then 
there is no swirl relative to the rotating frame observer. Therefore the 
inviscid cross force must be zero when q=v./r i# or equivalently when 

W=l ; This behavior is very similar to the quasi-static oscillation of an 
airfoil. The damping force there is related to the induced angle of attack 
due to the vertical motion. For the case of the whirling rotor, it is the 
induced inlet swirl angle change, due to the whirling motion, that 
creates the rotordynamic damping. Figure 6 shows the analogy between 
these two phenomena. 

This simple behavior for both the airfoil and the rotor is limited 

to the quasi static case. For the rotor, this condition of low reduced 
frequency, is satisfied when the fluid residence time in the seal is short 
when compared to the period of oscillation of the rotor. This can be 

expressed as 


7T <<1 ( 28 ) 

v x 

where V x is the axial velocity through the seal. An interesting and 
useful consequence of this relationship between the static and dynamic 
forces, is that in the absence of viscosity the direct damping, C„, can be 
calculated from purely static measurements of the cross stiffness , K„ y , 
versus the inlet swirl by 


C. =R, 


dK, y (V,) 

d(V.) 


( 29 ) 
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C„=R.i^ 


(30) 


coefficient fromTh? 1 bC , f USed 10 extract thc direct damping 

viscositv Hnwvr h h ° cocff,c ient measurements in the presence of 

viscosity. However, a discussion of this will be reserved until later 

The viscous contribution, F r , which is proportional to a'r and 
hence the swirl velocity difference Vj-V. contributes to the cross 

wSfn! ?" y ‘ ThC P hys ' ca] reason thal this force is independent of the 
" h ' r i" g i fre ? U f Cy ' S that 11 depends on a velocity difference which is 

whirl The mde «f endent K ° f any chan S e of reference frame due to the 
1. The mechanism that produces this force is shown Figure 7 The 

2 Sea ' Wilh a hi ‘ her s » irt «'oci,y ,ha„ Ss js: 
.., s ® a1, L f ss cxcess momentum enters through the narrow san on ihe 
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Dumn f U The nS ' d t C tHC ? land increasing thc static pressure like an ejector 

P- , P n Th p acc Wlth ,hc h'ghcst pressure will be at the bottom 

yielding a positive cross force as shown. A simple analysis of this is 
c y using the following simplified momentum equation 


^ (V " Vl) + ff = 0 


(31) 


this expression can be integrated to yield thc viscous cross force 

F _ ittq*R?(Vi -V*)r 

P h8* 
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If the swirl speed decreases through the seal, as is typical, then the 

frequency at which the force becomes negative shifts to a higher W, as 
shown. 

4. EXPERIMENTAL APPARATUS 

The Labyrinth Seal Test Facility (LSTF) was designed and built to 
measure the dynamic forces in a spinning/whirling labyrinth seal. 
Figure 9 shows a cross section of the hardware. Air from a compressor 
enters the first plenum and is turned radially outward through eight 1 
and 1/2 inch holes having honeycomb plugs. Next the air turns axially 
and accelerates through a set of replaceable swirl vanes into the swirl 
plenum. The air flows through the test seals and discharges to the 
atmosphere. The spinning/whirling motion of the seal is produced by 
the nested bearing arrangement. Different whirl amplitudes are 

obtained by adjusting the inner spindle bearing scat eccentricity. The 
spin motion (+6700rpm to -6700 rpm) is driven by an in-line flexible 
coupling driven by an electric motor. The speed of the whirling motion, 
which can be controlled independently, (+3400 rpm to -3400 rpm) is 
driven by a V-belt attached to another motor. Four equally spaced, flush 
mounted, high response, Kulitc XCS190 differential pressure transducers 
were placed on the seal land to measure the time resolved gland 

pressure oscillations created by the seal whirling motion. The back 
pressure ports were referenced to the gland average pressure to obtain 
higher sensitivity. Proximcters were used to precisely measure the 
whirling motion. Measurements of the swirl angle leaving the second 
knife were made with a hot wire anemometer. The data acquisition 
system was triggered and clocked with a chopper wheel attached to the 
whirl producing spindle. Thirty-two phase locked points were taken for 
each whirl revolution. Records of 64 revolutions were taken. The 
pressures were composite phase-locked ensemble averaged to find the 
forces acting on the rotor for each operating condition. A total of five 

different seal geometries, shown in Table 1, were tested under various 
operating conditions. The inlet pressure, swirl vane angle, spin 
velocity, whirl eccentricity and frequency were varied parametrically. 
Figure 10 shows a typical output trace from one of the Kulites versus 
data point number. Note the periodicity of 32 data points per cycle as 
expected. Spectral analysis of these signals showed that virtually all of 
the energy of the signal is concentrated at the whirl frequency for 
cases with no seal rotation [24], For cases including rotation, a small 
harmonic component, with amplitude of about 5% of the primary 
component due to whirl, at the spin frequency was present. This was 
due to very small deviations from circularity of the seal. 


S. COMPARISON OF EXPERIMENTS TO THEORY 

A strict comparison of the preceding theory with coupling to the 

LSFT test results is not possible because of the lack of control over the 

axial leakage gap, 5 C , the importance of which was shown in Section 2. 

What can be shown, however, is that for the range of likely values of 8 C , 

the upstream coupling effects increase the cross forces the required 
amount to explain the very large deviations between experimental data 
and the simple, uncoupled theory. 
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The general character of the comparison can be seen from 
Figure 11, which shows the nondimcnsional cross stiffness versus the 
nondimcnsional inlet swirl. The trends and comparisons for C„ arc the 
same as for K xy , due to the relationships shown in Section 3. The 
theoretical lines were calculated using the cavity parameters for Build 
I he volume ratio (upstream to seal cavity) was in that case 3.01 
whereas it was 4.64 for Build #2. The amplification factors calculated 

from the coupled theory (for 5 C = 10 mil) are 3.7 for Build #3 and 2 8 for 
Build #2. 

Looking first at the Build #3 data only, we notice a general 

grouping about the theoretical line for 8 C between 5 and 10 mil. The 
design value of this gap was 10 mil. We also see that the points appear 
in groups, each of which is aligned with a somewhat different 5 C line. 
These groups' in fact correspond to the different swirl vane assemblies, 
and o c can be expected to have remained constant within each of them, 
but perhaps not with exactly the same value from assembly to assembly. 

The data for the two outer groups of points for Build #2 arc 
roughly on a line with a slope lower by 1/1.3 than those for Build #3, in 
accordance with the noted difference between the volume ratios of both 
cases. The group for the smaller dimensionless inlet swirl is 

anomalously high, however. In general, for all builds, the data are 
indeed bracketed between the limits of the uncoupled (8 c -»°o) theory 
and the fully coupled (8 C = 0) models. 


As a related test of the theory, equations were derived for multi- 

cavity seals with upstream coupling [9]. This coupled theory for a 2- 

gland seal, plus the upstream cavity was compared to some 2-gland static 
sea data of Benckert and Wachter [5], whose facility had the same 
general upstream configuration as the LSTF. These researchers also 
ai ed to keep a careful control of their face seal, for which no value is 
reported. As noted before, no theoretical calculation had been able to 
match the higher than expected cross-forces in these short seals. Using 
he coupled model, the cross forces obtained in all three 2 gland seals 
that Benckert tested can be matched by the theory with appropriate 
choices of the axial scaling gap, (8 C =0) (a slightly different one for each 
build) Not only can the total force be predicted, but the relative 

contributions from each gland can be matched for the one ease for 

which this was reported. The model predicts that the effects of the 
upstream coupling die exponentially. After a few glands the upstream 

influence cannot be seen. The data of Benckert support this conclusion. 

Ih(> , 2 T1 , glVCS lhc dala of Benckerl a,on g with several eases from 

h n e ., model : The cross force as predicted with uniform upstream 
conditions is given. Finally, the value of the axial gap that matches 

sh™ that fh? , 1S 8iVCn f ° r CaCh CaSC - TheSC va,ues are reas °nable and 
show that the strong upstream coupling is the most likely cause for the 

high cross forces generated in the first gland of a seal. 

freauenfv 8 ^ f 2 3 lypical plot of the cross force vs. whirl 

lrequency for five different pressure ratios from Build #4 There is no 

rotation and the nominal inlet swirl angle is 8.6 degrees The general 

k dcpcndcnl behavior of the experiment da. i/ w 

Hnea^funton of T' y "?\ ^ parlicular * thc «*« force is a 
changes ^n m whir ' spccd K a "d the frequency at which thc force 
changes sign matches the theoretically predicted value within 
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experimental uncertainty. Of course the absolute force levels are much 
higher as predicted by the coupled theory. Figure 13 shows the cross 

force vs. whirl frequency for five different spin rates for the same 

build. It can clearly be seen that the data at different spin rates fall on 
parallel lines. Similar results were found for all builds at all flow 

conditions. The narrow implication of this is that changes in rotational 

speed do not affect the damping. More generally, these data strongly 

support the general validity of the cross force decomposition into ideal 

and viscous contributions as presented above. 

The method presented above for predicting the dynamic 
coefficients from inviscid static values can be readily extended to cases 
with "viscosity" if the frictional component is isolated properly. 
Probably the best method for removing this frictional "contamination" 
is to handle it directly in the experiment. If the inlet and exit swirl 

velocities are measured then for every value of the inlet swirl the seal 

spin speed should be adjusted to maintain a constant swirl difference, 
Vj-V*. Hence, the additive frictional component is fixed. For Vj = V\ the 

viscous component is totally eliminated. For typical designs the spin 
speed needs to be maintained at about 125% of the inlet swirl velocity to 
keep r = 0. If this procedure is not followed the measurements of a 
statically offset seal can still be used to obtain damping data. If static 
data is taken with no rotation, a correlation of the cross stiffness 
coefficient, K* y , versus swirl parameter, a, can be obtained. The 
problem is that as the inlet swirl velocity is increased the total cross 

force increases due to a higher inlet swirl (inviscid contribution) as 
well as the frictionally induced change in swirl. The change in the ideal 

component with inlet swirl, which is equal to the direct damping must 

be separated. One procedure is to use the static correlation to calculate a 
total force, F r (O d R,), at an inlet swirl velocity corresponding to a whirl 

frequency Q d . The viscous force, found from theory or correlation, is 
subtracted form the total force to yield the ideal component. The 

damping is the rate of change of this force with frequency. 

r F T (Q d R,)-F r (n d R,) F„(Q d R t ) 

The values of the direct damping, C^, calculated from the static 

measurement correlations employing this method agree with the 
values that were directly measured in the dynamic mode to within 13%, 
8% and 15% for builds 2,3 and 4 respectively as shown in Table 3. This 
demonstrates that it is possible to predict the dynamic coefficients with 
the use of static measurements only. Another method for extracting the 
dynamic coefficients from static data, based on extrapolation, is 
presented by Millsaps and Marlincz-Sanchez [26]. 

6. CONCLUSIONS 

A new theory was developed that contains the effect of flow 
coupling. It was shown that this theory is capable of predicting the 
larger than expected forces found experimentally. 

This finding has several important implications for design and 
analysis of short labyrinth seals, and, by extension, probably for other 
seal types as well. Data on cross forces from such seals should be 
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supplemented by a description of the upstream (and also the 

downstream) configuration of the test device, unless it can be 
ascertained by auxiliary tests that no coupling exists. As an example of 
the latter situation, it was verified that the rolordynamic force data 
were insensitive to the downstream configuration. This was found by 
blocking parts of the exhaust holes from the downstream chamber. The 
remaining holes were large enough to equilibrate this chamber with 
the atmosphere, hence making coupling negligible. 

Calculation of seal rotordynamic coefficients, either for design or 
for rotordynamic diagnostic purposes, should always account for these 
possible coupling to the external flowfield. Ignoring them may make 

other refinements, such as more precise 2 or 3-D cavity flow modeling 
or CFD calculations irrelevant by comparison. Including the upstream 
and downstream coupling is relatively straightforward for simple 
geometries, such as the LSTF rig, but no theory currently exists that can 
be used to predict the non-uniformities ahead of, for example, the seal 
on a turbine lip shroud, where the tip area interacts strongly with the 
main turbine through flow. Development of such a theory would seem 
important. Similar considerations should apply for other seal 

environments. 

In future designs of rotordynamic testing devices, consideration 
should be given to cither minimizing the coupling mechanisms or 
matching the expected coupling levels that may occur in a real 

lurbomachinc. In any case measurements of the flow field upstream and 
downstream should be done to asses the degree of non-uniformity 
induced from the offset seal. 

The following conclusions on the nature of rotordynamic 

damping have been drawn from consideration of the analytical model 
and the supporting experimental data. 

1. The total cross force acting on a labyrinth seal at a given whirl 

frequency can be decomposed into "ideal" and "viscous" components. 

2. The ideal component, which is due to an inviscid swirling flow, is a 

unique function of the inlet swirl relative to the gap variation phase 

speed. This force component vanishes when the velocity of the 

traveling gap is equal to the swirl velocity inside the gland. This 

component is solely responsible for damping. 

3. The viscous component docs not create nor alter the direct damping. 
It adds to (if the swirl velocity decreases in the gland) or subtracts from 
(if the swirl increases) the cross stiffness only. 

4. The direct damping can be calculated from measurements of cross 

stiffness. The importance of this is that difficult and expensive dynamic 

measurements arc not necessary in order to obtain damping 

coefficients. 
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Figure 1. Schematic of Labyrinth Seal Test Facility test section showing 
upstream swirl cavity and the flows into and out of the center hub 
plenum. 



Figure 2. Ratio of rotordynamic damping with upstream coupling to 
that with no coupling versus swirl cavity to seal gland area ratio. 
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Figure 3. Predicted direct and cross force at £2 = 0 for the coupled model 
versus the relative axial clearance. The geometry is the same as for 
build #3. Jl s =1.4, a v = 15 degrees, Ej =0.1407 and 0)=0. 



Figure 4. Direct and cross force predicted by the coupled model versus 

the center cavity swirl velocity. These are for the same conditions as 
Figure 3. 
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. Effect of relative whirl amplitude on the direct and cross 
as predicted by the coupled model with hub leakage. Build #3 

and 8 c /8j=0.19. 
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airfoil 


6 Analogy between quasi-static damping of a translational 
generating an induced angle of attack and a whirling seal 


inducing a inlet swirl component. 
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7 * ,f hysical origin of the v ‘scous force. This force is due to a 

ThAr C r y d, . ffe ^ encc and hence 'S independent of the whirl frequency 
Therefore, it does not contribute to the damping. 4 * 
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" General frequency dependent behavior of the 
both the ideal and viscous contributions. 
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Figure 9. Cross section of the LSTF test section and nested spindle 
whirl/spin producing rotating machinery. 



Figure 10. Typical output trace from one of the flush mount pressure 
transducers. There are 32 data points per whirl revolution. 
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THEORY EXPERIMENT 



Figure 11. Nondimcnsional cross stiffness, K xy , versus the swirl 

parameter, a, for the experimental data and theory, with the axial gap 
used as a parameter. All experimental values fall between the 
theoretical predictions with 0.004" (0.0001m)<8 c <0.017" (0.0004m). The 

top thick line is the fully coupled case (i.e. no center leakage) and the 
bottom on is for uniform inlet conditions (i.e. no coupling). Calculations 
are for build #3 geometry. 



Figure 12. Experimentally obtained cross force versus whirl frequency 
tor five different pressure ratios. These data are from build#4 with 8.6 
degrees of inlet swirl and no shaft rotation. 
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Table 2. Comparison of the static data of and Wachter to the coupled 
model. The first column shows the experimental value. The second 
column gives the value predicted with full coupling. The third gives 
predictions for constant upstream conditions (i.e. no coupling). The last 
column gives the value of the axial sealing gap needed for the model to 
match the experimentally obtained value. 


CONFIG. 

Fr(meaj,) 

MK = 0) 

Ft{S* = oo) 


1 

10.21(N) 

16.28(N) 

4.99(N) 

0.008" 

2 

8.28(N) 

16.15(N) 

4.25(N) 

■ 

3 

ll.91(N) 

15.51(N) 

4.09(N) 

■ 


Table 3. The first column gives the cross stiffness correlation 
coefficient obtained from static measurements. The next two columns 
give the total, and frictional force, respectively calculated at a whirl 
speed of 300 (rads/sec). The last two columns give the calculated and 
measured direct damping coefficients. 


BUILD# 

dc 

F T [O t ) 

H 

B 

From Static 
Correlation 


2 

0.372 

7763 

953 

0.289 

22.70 

19.95 

3 

0.416 

27821 

3162 

D 

82.19 

75.80 

4 

0.283 

20001 

3162 

0.247 

56.13 

48.54 

5 

0.338 

7053 

953 

0.231 

20.33 

15.81 
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